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Abstract. This article was written to graduate and postgraduate students 
of physics. We intend to be didactical using a mathematical approach in 
Riemannian spaces as rigorous as possible. Brief comments are done on 
topics of General Relativity (GR) and Special Relativity (SR). 
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(I) Introduction. 

Are essentially shown mathematical aspects on Special Relativity, 
Riemannian Geometry and General Relativity. Comparison with 
experimental results can be seen in references mentioned in the text. In 
Section 1 we see Manifolds: coordinate transformations, metric and affine 
spaces. In Section 2 are given examples of metric spaces and their 
geometries. In Section 3 we define Geodesics. In Section4 are analyzed 
local Cartesian or local inertial coordinate-systems. In Section 5 are 
shown effects due to the Earth gravitational field (Schwarzschild Metric). 
In Section 6 we see the Equivalence and the Geodesic Principles. In 
Section 7 we consider accelerated coordinate-systems in SR. In Section 8 
we calculate spatial distances and time intervals in GR. 
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(1) Manifolds. 

In elementary analytical geometry 1131 the fundamental notion is that 
geometrical points are defined by means of its coordinates. Thus, in the 
geometry of the Euclidean plane, for instance, a point is specified by 
giving its two Cartesian coordinates (X,Y) or its polar coordinates (r,0). 

All points which together constitute the plane are said to form a 2- 
dimensional manifold or 2-dimensional space of points, the number of 
dimensions of the manifold being equal to the number of independent 
coordinates required to specify a point in it. Generalizing this idea, an n- 
dimensional manifold or n-dimensional space of points is one for which n 
independent real numbers (x^x 2 ,...^ 11 ) are required to specify every point 
completely. These n numbers are denoted collectively by {x 1 }^ and are 
called the coordinates of the point. For the moment the manifolds are 
assumed to have no structure except that it is continuous in the sense that, 
in the neighborhood of every point (x) there are other points whose 
coordinates differ infinitesimally from those of (x). Such neighboring point 
has coordinates (x + dx), the small quantities (dx 1 , dx 2 ,...,dx n ) being called 
the differentials of the coordinates (x). Being these coordinates 
independent we have dxVchC = 5/. 

(l.l)Coordinate transformations. 

Let us consider two manifolds X ={x 1 } i=L .. n and X'={x' J }j =L .. m .To 
perform the coordinates transformation X —> X' it is necessary that the 
manifolds have the same dimensions, that is, n = m. The operation by 
which the coordinates (x) of every point in the manifold are altered to (x") 
is called coordinate transformation. The coordinate transformation X —> X' 
are defined by n equations x A = f '“(x 1 , x 2 „...,x n ) where X = 1,2,3,...,n and 
the functions f x are soluble in such a way that would be possible to obtain 
the inverse transformation, X'— ► X, that is, x x = g > “(x' 1 , x' 2 ,,...,x' n ). [1 ‘ 3] 
Differentiating partially the x A (dx v ) we obtain the transformation law for 
the differentials 

dx A = (df Vdx v )dx v = (dx'Vdx v )dx v (1.1.1). 

Mathematical conditions that this transformation be (theoretically) possible 
is that the determinant of the coefficients ( Jacobian ) must not vanish 
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detlldx'VdxJI f 0 


(1.1.2) 


at all points, or nearly all points. If this inequality does not hold, then the 
equations X'={x 4 }j = i ... m do not a legitimate this set of coordinates. [4] 

(1.2) Space Affine: Co variant and Contravariant Tensors and Scalars. 

Any aggregate of n quantities IF which transform from S to S' like 

U' v = (dx' v /&giF (1.2.3) 

is called a contravariant vector ( 1st order tensor ). Quantities which 
transform like 

V' v = (0x70x' v )U v (1.2.4) 

are called covariant vector. Although covariant and contravariant vectors 
have nothing to do with each other their inner product defined as IF is 
an invariant, that is, it is independent of the coordinate system. Indeed, 

V V U 41 = (0x70x'Q (0x'70x p )V a U p = (0x a /0x p ) V a U p = 5 P 8 V a U p , 

that is, 

VV U 41 = V a U a = Invariant (1.2.5), 

showing that their inner product is a scalar. 

Consider a function F(xQ = F(x' v ). Since dF(xQ = dF(x' v ) the 
differential dF is an invariant, that is, it is a scalar. 

A space in which covariant and contravariant vectors exist separately 
is called affine. In other words, in an affine space a vector does not exist; 
but covariant and contravariant vectors exist independently of each other 
Such spaces are very general and have been extensively studied by 
mathematicians. [5] The scalar is Oth order tensor and U v or U v ar e first 
order tensors. More details about definitions and properties of covariant 
and contravariant tensors of any order can be seen elsewhere. [151 For 
instance, a quantity U pv behaves like a 2. nd order covariant tensor if it 
transforms as 

UV = (dx a /dx'Q (0x p /0x' v ) U a(3 (1.2.6). 

Higher order tensor can be seen, for instance, in references 1-5. In 
GR there are two important tensor, one of rank four R\ v ^ called Riemann- 



Christoff el tensor and one of rank two named Ricci tensor. These are 

shown explicitly in Appendix A. 

(1.3) Metric Spaces. 

For the moment the manifolds were assumed to have no structure 
except that it is continuous in the sense that, in the neighborhood of every 
point (x) there are other points whose coordinates differ infinitesimally 
from those of (x). In particular no definition of distance ds between pair of 
points has been given. In arriving at such a definition is introduced a metric 
that contains within itself the essentials of the geometry of the space. 

In fundamental physical theories we deal with objects or quantities 
which are independent of the particular choice of the mode of description. 
In physics it is desirable to deal with a space in which the concept of 
distance and tensors exist as an objective physical reality. Such spaces are 
called metric. Before to define rigorously metric spaces let us remember 
that in physics the basic observable geometric quantity is the distance ds 
between two points. For instance, the most familiar metric space is 3- 
dimensional Euclidean space. In fact, a "metric" is the generalization of the 
Euclidean metric arising from the four long-known properties of the 
Euclidean distance. In a 3-dim Euclidean space, in Cartesian coordinates, 
(x 1 , x 2 , x 3 ) = (x, y, z) the infinitesimal distance ds between two points is 
defined by the invariant ds 2 = dx^dx^, that is 

ds 2 = dx 2 + dy 2 + dz 2 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 = gy dx‘dx j , (1.3.1) 

where g 1J = 8y is called " metric tensor”. The distance ds between the same 
points must equal to that given by (1.3.1) independently of the coordinate 
system adopted. For instance, in spherical coordinates x' 1 = r , x' 2 = 9 and 
x' 3 = (p we must have [1 ' 5] 

ds 2 = ds" 2 = r 2 + r 2 d0 2 + r 2 sin 2 0 dcp 2 = g'y dxAbC j (1.3.2), 

2 2 2 

where now the metric tensor is given by g' l} = (1, r , Csin“0). 

Note that to get the above results the transformation X —> X' is given 
by x = r sin0 coscp, x“ = r sin9 sincp, x = r cos0. The distance ds 
between two points must be independent of the coordinates system. 
Similarly, vectors like velocity V, force F, distance r, ...must be 
independent of the coordinates used to represent them. Let us see how 
these ideas can be generalized involving n-dimensional manifolds. Before 
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to do this note that ds 2 in the 4-dim in the SR pseudo- Euclidean 
Minkowski space x 3 = (x',x 2 ,x 3 ,x 4 = x°= ct) is written in different ways in 
text books and papers. That is, we have ds 2 = (x 1 )+(x 2 ) 2 +(x 3 ) 2 - c 2 dt 2 
where gn = g 2 2 = g 33 = 1, gik = 0 (#k) and g 0G = -1 or - ds 2 and, sometimes, 
ds 2 = c 2 dt 2 -{(xV(x 2 ) 2 + (x 3 ) 2 } with g n = g 22 = g 33 = -1 g ik = 0 (#k) and 
goo = 1. Since ds 2 is invariant this equation is valid in all reference system. 
An interesting case is a clock at rest in a reference S 0 .Since it is at rest its 
displacements vanish, dx 0 = dy G = dz 0 = 0. Let us denote the time measured 
by this clock by dt = dr, where the time x measured in S G is called proper 
time, so ds~= - c dx~. Now, let us consider a system S which moves with 
constant velocity v relatively to S Q . An observer in S will observe that this 
clock in a time interval dt suffer displacements dx, dy and dz in S. So, we 
can put ds 2 = - c 2 dx 2 = dx 2 + dy 2 + dz 2 - c 2 dt 2 . In this way, taking into 
account that v 2 = (dx 2 + dy 2 + dz 2 )/dt 2 , we obtain dx = [1 - (v/c) 2 ] 172 dt, as 
seen in basic SR courses. [6] 

There are spaces in which covariant and contravariant vectors do not 
exist independently but they can be converted into each other. These have 
a further property that the index of a contravariant vector lA can be 
lowered to become by the operation 

u,= g, v U v (1.3.3). 

The reverse transformation is would be given by 

U v = g v ^U, (1.3.4), 

thus, U, = g^ v U v = g^ v g va U a which means that g^ v g va = 8^ a . In other 
words, g^ v is the inverse of g^ v , and vice-versa, that is, g^ v = M^/lgl, where 
Igl is the determinant of g^ v and is the minor of the element g^ v . In a 
metric space a vector is represented by U and has an objective meaning, 
independently of the adopted mathematical representation. Identical 
properties are obeyed for higher order tensors,covariant or contravariant 15] 
The square of a vector U" is an invariant defined by the inner product 

U 2 = = g, v lTTJ v = g aP U a Up (1.3.5). 

Similarly, we define the square ds 2 = ds 2 of the differential dx 11 by: 

ds 2 = ds 2 = g^ v dxAlx v (1.3.6). 
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The quantity ds (1.3.6) is called the line element of the metric space. It 
expresses the invariance of ds , also named fundamental quadratic form of 
metric spaces. An important assumption which is not explicitly expressed 
in usual textbooks, is that the basic observable of the metric geometry is the 
distance ds. l4] Among other things, we can deduce from (1.3.5) that g^ v is a 
symmetric tensor, that is, g^ v = g Vi ,. For instance, Riemannian spaces used 
in General Relativity (GR) are metric. In most of the applications of these 
Riemannian spaces, special types of spaces (said to be diagonal) are 
considered in which g^ v = 0 if p f v. 

It can be proved that the functions g^ v are components of a covariant 
tensor of rank two. Indeed if in the coordinate transformation X —> X': 

ds 2 = g^ dx A dx^ = g ^ dx\bC , 
using (1.2.6), that is, taking dx A = (dx'VSx^dx^ we get 

g 'a.p (xo = g „(x')(axVax' a )(ax75x' p ), ( 1 . 3 . 7 ), 

where x 1 = x 1 ({x}) showing that g^ v transforms as a 2. nd order covariant 
tensor according to (1.2.6). Vice-versa we must have when X'—► X : 

gap = gnv({x J }) (dxVdx a ) (dx' v /dx p ) (1.3.8). 

Remembering that the space-time metric is defined by the line 
invariants^ 3] ds 2 = g^ v dx 11 dx v , g^ v e x^ are, respectively, the metric tensor 
and the coordinates of the referential system. 

Note that in all general "real" space-times the determinant Igl of the 
gi k must be Igl < 0. [1,3] The gravitational field is named constant when all 
components of the tensor g a p do not depend of the temporal coordinate x 0 . 
When the components g 0 p = 0 the champ is said to be static. It is stationary 
when the components g 0 p f 0. 

(2) Examples of Metric Spaces and Geometries. 

Fundamental aspects that are very important to be mentioned are 
those given by relations between the line elements ds. By coordinate 
transformations different geometries be generated. [2 ’ 3] 

(2.1) Two and Three-dim Euclidean spaces. [3] 

In a 2-dim Cartesian orthogonal coordinates (X = x\Y= x 2 ) the 
metric ds is given by, ds 2 = dX 2 + dY 2 (2.1.1). 
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This formula, where g 11 = g 22 = 1, that expresses the metric of the space 
defined by the variables (X,Y) contains the essentials of the Plane 
Euclidean Geometry, straight line is the shortest distance between two 
points, parallel lines are of infinite length and do not intersect in any finite 
part of the plane, and so on... These properties that are due to the use of 
Cartesian coordinates disappear if other types of coordinate-systems are 
employed. For example, if we transform from the Cartesian coordinates 
X = x and Y= x~ to polar coordinates r = x' , 0 = x' by the equations 
x^x^cosx' 2 and x 2 = x' 1 sin x' 2 the metric (2.1.1) becomes 

ds 2 = dr 2 + r 2 d0 2 (2.1.2), 

with g 11 = 1 e g 22 = r 2 . It does not describe a Plane Geometry. 

In 3-dim Euclidean space, when Cartesian coordinates (X,Y,Z) are 
used, the metric has the form 

ds 2 = dX 2 + dY 2 + dZ 2 (2.1.3), 


which is again the statement of Pythagora's theorem. If these coordinates 

12 3 12 3 

are changed as (x = X, x~ = Y, x = Z) —> (x' = r, x' = 0, x' = cp) according 


. I A ■ A A 

to x = x sin x cos x , 


2 A ■ A ■ A i 3 A A 

x = x sin x sin x and x = x cos x we 


see that (2.1,3) becomes written as 


ds 2 = dr 2 + r 2 d0 2 + r 2 sin 2 0 dip 2 (2.1.4). 

2 2 2 

Taking r = constant = a the metric becomes (with gu = a and g 2 2 = a“sin“9) 

ds 2 = a 2 (d0 2 + sin 2 0 dcp 2 ) (2.1.5), 


that describes the geometry of a 2-dim surface of a sphere with radius r = a 
which is intrinsically different from that of the 2-dim Euclidean plane: 
straight lines are replaced by great circles with finite length, and 
necessarily intersect, with no parallels in the Euclidean sense, and so on... 

Note that these two sets of equations in 2 and 3 dimensions have the 
property of being soluble for the (x") in terms of the (x); for example, 2- 
dim case the equations we have the relations 

x' 1 = {(x 1 ) 2 + (x 2 ) 2 } 172 and x' 2 = tan' 1 (xV). 


In the 3-dim case we have similar relations, but more complicated. 
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(2.2) 2-dim Coordinate-System (u,v) where X = uv and Y = (u 2 +v 2 )/2. 

In this case ds 2 = dX 2 + dY 2 becomes written as 

ds 2 = (u 2 + v 2 )du 2 + 4uvdudv + (u 2 +v 2 )dv 2 , (2.2.1), 

with gn = u 2 + v 2 , gi 2 = g 2 i = 4u and g 2 2 = gn= u 2 + v 2 , respectively, where 
the orthogonality property is masked in the new coordinate-system. What 
kind of geometry does this (u,v) coordinates describe? This question is very 
difficult to be answered in a general case. 

(2.3) 3-dim manifold (X,Y,Z) with the metric 

ds 2 = dX 2 + dY 2 - dZ 2 (2.3.1), 

(gn = g 2 2 = 1 and g 33 = -1) is a 3-dim Minkowski space different from the 
Euclidean one. If the coordinates are transformed to (ri,0,tp ) like 

X = r| sinhG coscp , Y = ri sinhB sincp e Z = r! coshO 

we verify that (2.3.1) becomes 

ds 2 = - d ri 2 + fl 2 d9 2 + fl 2 sinh 2 0 dcp 2 (2.3.2). 

Taking r = a = constant we obtain ds 2 = a 2 (d0 2 + r^ sinh 2 0 d(p 2 ). With a 
new coordinate transformation 0 —> r, defined by sinh0 = r/(l- r /4) we get 

ds 2 = a 2 (dr 2 + r 2 dq) 2 )/(l- r 2 /4) 2 (2.3.3). 

From (2.3.3), r = 0 when 0 = 0 and r = 2 when 0 —> oo. Again, it is difficult 
to determine what kind of geometry is defined by (2.3.1)... 


(3)Geodesics. [3] 

In preceding sections have been performed changes of tensors when 
the coordinates of the points of a Riemannian space are transformed. Now 
we consider changes of a different kind that arise when the coordinate - 
system being kept fixed and the value of the tensor at one point is 
compared with its value at another point. Such changes may be usefully 
regarded as due to the "motion" of the tensor from one point to the other 
and it is therefore necessary to define some kind of path through the 
Riemannian space along which the tensor may be "imagined" to travel. 
These fundamental paths are called geodesics of the space and they have 
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properties analogous to those of straight lines in Euclidean space. The 
geodesics are particular kinds of curves in the space defined by n equations 

x l = F>), (I=l,2,3,...,n) (3.1), 

where u is parameter varying from point to point of the curve. Substituting 
equations (3.1) into (1.3.6) ds 2 = g^ v dxMx v and then integrating with 
respect to u, it is possible to express the interval s measured along the curve 
in terms of u, that is, 


S= {g (x)(dx /du)(dx / du)} 1/2 du 

0 (3.2), 

where u G and Ui are the values of p at the points P Q and P 1? respectively. The 
geodesic joining the two points P Q and Pl is then defined to be a curve for 
which the interval s between P Q and Pj has a stationary value [4] ("stationary 
interval ") that is, a 4-dimensional "length". Stationary here is used in the 
sense in which that term is used in the calculus of variations, namely, that 
the interval s measured varies minimally along any other neighboring curve 
joining the two points. Any other curve joining P 0 and , and always lying 
close to the geodesics, will have equations of the form 

x*'“ = x x + 8 C 0 X = F'“(u) + so/fu) (3.3), 

A, 

where co = 0 at u = u G and u = ui, and s is a small quantity whose square 
and higher powers may be neglected. The new interval s* calculated along 
the neighboring curve joining P 0 and Pi would be given by 

i 

S = {g (x )(dx /du)(dx /du)} 1/2 du 

o (3.4). 

Calculating the difference s*- s neglecting all powers of g above the first 
order one can verify that [3] to this interval to have a stationary value for the 
geodesic compared with the neighboring curves, s*- s must be zero for any 
choice of the functions co . In this condition one can prove that the 
following condition must be satisfied: 

d 2 x T /ds 2 + T^ V T (dxVds) (dx v /ds) = 0 (x = l,2,3,...,n) (3.5), 
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where I\ V T = (kp,v) = (pk,v) = (1/2) { cg-Jcsd + dgjdx x - dg h Jdx v } 
are the Christoffel symbols [3,5] . Equations (3.5) are the standard forms of 
the geodesics equations of the Riemannian space, excluding, however, 
those geodesics for which s = 0 along the curve. As the manifold has n 
dimension, the geodesic equations are a system of ordinary differential 
equations for the coordinate variables. If dx M correspond to an infinitesimal 
displacement along the geodesic for a change ds of interval, the vector n 11 = 
dx'Vds is called the unit tangent vector to the geodesic. 

By dividing (1.3.6) by ds 2 it follows that 

g^vidx'VdsXdxVds) = 1 (3.6), 

which shows that n 11 is a unit vector (vector of unit length). Another 
important conclusion [3] is that (3.6) is an integral of the n equations of the 
geodesic (3.5)(see also Appendix B). 

According to Appendix C a vector q^ is timelike, null or spacelike: 

timelike when g^ v (x) q^ q v > 0, 

null when g ? , v (x)q >v q v =0, (3.7) 

spacelike when g ? , v (x)q >v q v <0. 

(3.1)Timelike geodesics. 

The equation of motion of a particle with mass in GR is written as 
Dp'Vds = T 1 , where D is the covariant derivative (see Appendix C), 

p 11 = mv^ , v^ = dx'Vds, ds 2 = g^ v dx |A dx v and F is the force acting on the 
particle, that is, 

dWds 2 + r,/ (dxVds) (dx v /ds) = f (3.1.1). 


In the case of free particle for which f 11 = 0 , Eq.(3.1.1) reduces to 


d V/ds 2 + r } / (dxVds) (dx v /ds) = 0 (3.1.2), 


which reinforces that the path of a free massive particle is a geodesic 
according to Eq.(3.5) in the spacetime. This result is often stated as a 
explicit postulate of GR, known as geodesic postulate, but emerges here as 
a natural consequence of the way in which we generalize SR concepts. [8] 
The path of a free massive particle in flat spacetime of the SR is a straight 
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line and this generalizes to a geodesic in curved spacetime of the GR. 

If, besides (3.1.2), is obeyed the condition g^ (dx |i /ds)(dx7ds) = 1 we 
say that the geodesics is timelike. At any point on the path of a massive 
particle its world velocity is a tangent vector (n^) to the path; Eq.(1.3.6) tell 
us that this tangent vector is timelike. So, this particle follows a timelike 
path through the spacetime, and in particular a free particle follows a 
timelike geodesic. 

(3.2) Null Geodesics. 

A null geodesic, is obtained by assuming that the interval ds = 0 
between any two points on the curve. For instance, from (3.4) we obtain 

g^ v (x)(dxVdu)(dx v /du) = 0 (3.2.1). 

This equation in SR generalizes the relation rj^ v (x)(dxVdt)(dx v /dt) = 0 , 
where r^ = diagf- 1,1,1,1] and the velocity v = (dx /dt). This is equivalent 
to c 2 dt 2 - dx 2 - dy 2 - dz 2 = 0. Trajectories of light signals in 3-dim Euclidean 
space are described by null geodesics that are straight lines. In this space 
ds 2 = dX 2 + dY 2 + dZ 2 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 = gy dx'dx', where the 
" metric tensor " g 1J = 5y. As g lj are constants from (3.5) we have d 2 x T /ds 2 = 0, 
that is, d 2 X/ds 2 = d 2 Y/ds 2 = d 2 Z/ds 2 = 0 which are immediately integrable to 
give the equations of the straight line, passing by the point (X 0 ,Y 0 , Z 0 ): 

( X-X 0 )/f = (Y-Y 0 )/m = (Z-Z 0 )/n = s, (3.2.2) 

where the direction-cosines dxYds = cosG^ (the unit tangent vectors) of the 
line are (cos9 x = i, cos9 y = m, cos9 z = n) which satisfy the condition, 
following (3.6), 

f 2 + m 2 + n 2 = 1. (3.2.3). 

In Riemannian space null geodesics, instead of (3.5), are given by 

g^ v (x)(dxVdu)(dx v /du) = 0 (3.2.4) 

and 

d 2 x T /du 2 + r^ v T (dxVdu) (dx v /du) = 9 (x = 1,2,3,...,n) (3.2.5). 

For a photon that follows a null geodesic Eq.(3.2.4) tell us that the tangent 
vectors to its path are null. In this case it is more adequate to write (3.2.5) 
using the fact that the direction of the light propagation in geometric optics 
is given by the wave quadri-vector k T = dx T /du tangent to the light ray. 

Thus, (3.2.5) becomes 

dkVdu + T^ V T k^ k v = 0 (x = l,2,3,...,n) (3.2.6). 
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(3.3)Spacelike Geodesics. 

Spacelike paths and spacelike geodesics may also be defined 
following the procedure adopted above, but these have no physical 
significance unless one believes in tachyons. [8] 

Another Examples of Geodesics 

(1) Null geodesics in a 3-dim Minkowski space (X,Y,Z) with the metric 
ds 2 = dX 2 + dY 2 - dZ 2 

In this case as gn = g 2 2 = 1 and g 33 = -1, following the procedure 
adopted in Section (3.2) we obtain 

i 2 + m 2 - n 2 = 0, instead of (3.2.3). 

(2) Great Circles on the Surface of a Sphere. 

In a sphere ds 2 is given by (2.1.4): ds 2 = dr 2 + r 2 d0 2 + r 2 sin 2 0 dcp 2 . 
Curves with cp = constant are the great circles ("meridians"); curves with 
0 = constant are the small circles of latitude; only the 0 = n/2 is a great 
circle. If in the geodesic equations d 2 x T /ds 2 + Tj w T (d 2 xVds)(d 2 x v /ds) = 0 we 
put x = 0, x = cp and r = a = constant these equations reduce to the pair 

d[a 2 (d0/ds)]/ds -a 2 sin0cos0(dcp/d0) 2 = 0 and dfa 2 sin 2 0(dcp/ds)]/ds = 0. 

Particular solutions of these equations are obviously: 

(a) cp = constant and a0 = s, corresponding to the "meridians". 

(b) 0 = 7i/2 and acp = s, corresponding to the "equator". 

Note that circle of latitude on which 0 = constant 7 tt/ 2 does not 
satisfy the geodesic equations. This implies that great circles are geodesics 
and small circles are not geodesics. 


(3.4)Lagrangian formalism and Geodesics. 

Defining the Lagragian function [8] L = L(x a , x a ), with x a = x a (u), 

L(x a , x a ) = (l/2)g, v (x a ) (dx*7ds) (dx v /ds) (3.4.1), 

where x a = (dx a /ds), the Euler-Lagrange equations are given by 

d(0L/ax7/du - oL/ox M = 0 (3.4.2). 

From (3.4.1) and (3.4.2) one can deduce the geodesic equations (3.1.2), 

dV/ds 2 + D/ (dxVds) (dx v /ds) = 0 . 

This implies the geodesics, timelike or null, can be directly obtained using 
(3.4.1) and (3.4.2). In Appendix D , using the Lagrangian formalism, we 
see how to calculate trajectories of particles and photons in the vicinity of 
spherical massive objects. 
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(4) Local Cartesian or Local Inertial Coordinate-System. 

In a Riemannian space of n dimensions there are certain types of 
coordinate-systems (X) which are useful for describing the points lying in 
the neighborhood of a given point 0. [3] One example is the Local Cartesian 
Coordinate-System (LCCS) or Local Inertial Coordinate-System (LICS) 
which exists at any point of a Riemannian space of any type. We consider 
here only the case when the metric is orthogonal, that is, 

ds 2 = £/Ju(x)(dx x ) 2 (4.1), 

where g u = £>.yu(x) and £>, are equal to either +1 or -1. Let the point O 
have coordinates x G and consider the coordinates 

xMyuWI'V-Xo 1 ) (X = 1,2,3,..., n) (4.2), 


where the summation convention being suspended. The differentials of 
these coordinates are dX A = [yu(x„)] l/2 dx A and therefore the substitution 
into (4.1) gives the metric 

ds” — I^dX 1 ) 2 (4.3), 

which is identical with the metric of a flat space. [3] The coordinates X ' 
are called local Cartesian and are valid for those points near enough to O 
for the differences [yu(x)] [yu(x () )] “ to be of order not exceeding 

(x^ - x Q x ). If the space were indeed flat, then it would be possible to 
integrate the differential equations dXVdx 1 = [yu(x)] 1/2 , a process that 
cannot be performed in a curved Riemannian space. 

As is shown in reference [!] in LICS the Christoffel symbols T^ 1 can 
be annulled, that is, Txv = 0 and, consequently, the geodesic equations (3.5) 
become d 2 x T /du 2 = 0 (x = 1,2,3,...,n). The LIC are also called locally 
geodesic coordinates. 

(4.1)Local Inertial Coordinates in Schwarzschild Space-Time. 

The metric in the polar Schwarzschild Spacetime (SST) is given 
by, [1 ' 3 ’ 5] where % = GM/c 2 , 

ds 2 = (l-2y/r)c 2 dt 2 - dr 2 /(l-2y/r) - r 2 d0 2 - r 2 sin 2 0d 2 (p (4.1.1). 

Putting r = (1 + j/ 2r*)r*, ct = x 4 , x J = r*sin0, x 2 = r*sin0sin(p, x 3 = r*cos0 
the metric ds 2 given by (4.2), at the event (x^x^, x 0 2 ,x 0 3 ) introducing the 
LIC defined by 

X 4 = [(1 -x/2r 0 *)/'(1 +x/2r 0 *)] (x 4 -x 0 4 ) and X 1 = (l+x/2r 0 *) 2 (x' - x 0 ‘), 
becomes [3] 

ds 2 =(dX 4 ) 2 -Xi(dX 1 ) 2 (4.1.2). 
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In the case of the Minkowski Space -Time (MST) where the metric is 
given by ds 2 = (dx 4 ) 2 - E (dx 1 ) 2 the velocity qj of the particle is defined by 
qi = cE (dx 1 ) 2 ] 1/2 /dx 4 . [3] Hence at the event (x 0 ), the local Cartesian 
velocity q G of the particle may be defined in manner analogous by 

q„ ={c[I, (dX‘) 2 ] 1,2 /dX 4 ) 0 = [(l+x&o^^l-x&o^KScfdxf/dx 4 ),, (4.1.3). 

LIC may be set up at any event and one may also work back to the polar 
system (t,r,9,(p). 

(5) Earth Gravitational Field (Schwartzschiid Metric). 

Let us assume that the Earth is a spherical body with mass M, radius 
R and with angular velocity Cl. Taking the Earth center of mass as the 
origin of an inertial system S one can show that the gravitational field 
generated by the Earth in vacuum (disregarding the Earth's spin effect) is 
described by the Schwartzschiid metric(SM) [13 51 ds 2 that in spherical polar 
coordinates (r, 9, cp) is given by 

ds 2 = (1-2%/r) c 2 dt 2 - dr 2 / (1-2%/r) - r 2 d9 2 - r 2 sin 2 9d 2 (p (5.1), 

where % = GM/c 2 . It is always emphasized that in GR the variables t and r 
are coordinates, that is, they are not necessarily "physical quantities" such 
as time and distance in classical mechanics. In Appendix E we analyze the 
Length and Time in the SST. 

(5.1) Gravitational Doppler Effect. 

Consider two clocks 1 and 2 at rest with respect to Earth with 
coordinates r! and r 2 . As for both dr = dep = d9 = 0 using (5.1) we obtain 
ds = c“dx where we made t = x, defined as proper time 1 that is, the 
time measured in the reference frame in which the clock is at rest. Thus, for 
points 1 and 2,we get 

dii = Vg 00 (ri) dt = (l-2Gm/c 2 ri) 1/2 dt 

dx 2 = Vgoo (r 2 ) dt = (l-2Gm/c 2 r 2 ) 1/2 dt (5.1.1) 

When r —> oo the proper time x becomes equal to the time coordinate t. 
Using (5.1) we can compare the time intervals that would be measured in 
reference frames 1 and 2. According to (5.1.1), if r 2 > r t the clock in 2 
moves faster than the clock in 2. Time measurements made with atomic 
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clocks of high precision have confirmed the predictions made by (5.1.1). [7] 
In these clocks the time measurements are based on the frequencies of 
vibrations v of atomic transitions. Thus, in 1 and 2 we would have, 
respectively, the frequencies V|=l/di| and v 2 =l/dx 2 given by, where v =l/dt, 

Vi = v/(l-2Gm/c 2 r!) 1/2 and v 2 = v/(l-2Gm/c 2 r 2 ) 1/2 

that is, 

vl/v2 = {(l-2Gm /c 2 r 2 ) / (l-2Gm /c 2 ^)} 172 (5.1.2). 

Equations (5.1.2) show the effect of the change of the 
electromagnetic frequencies generated by the gravitational field. It is the 
Gravitational Doppler Effect. 

(5.2) Gravitational & Transverse Kinematic Doppler Effect. 

To evaluate the change in frequency v due to the gravitational field 
and to the speed of movement of the clocks let us assume that one clock is 
fixed on Earth and the other is rotating around the Earth. Take clock 1 fixed 
on the surface of the Earth. Thus, 1 has a rotational motion with angular 
velocity Q relatively to the inertial system X. Let us assume that 2 is fixed 
on an airplane that carries out a circumnavigation movement around the 
Earth (for example, along the terrestrial equator) with a constant tangent 
velocity V. The fixed clock 1 on Earth (reference 1) putting Q = d9/dt and 
dr/dt = 0, using (5.1) would measure a time interval di = dx given by 

(dx T ) 2 = {(1-2GM/c 2 R) - R 2 Q 2 /c 2 }dt 2 (5.2.1). 

It is important to note that in (5.2.1) besides the gravitational effect 
given by the term (1-2GM /c“R)“dt“ there is also a kinematic effect given 
by - (R 2 G 2 /c 2 ) dt 2 which in SR is responsible for the Transverse Doppler 
Effect. For the fixed reference 2 in airplane we have r 2 = R + h, where h is 
the height of the flight; as the airplane moves with velocity V relative to the 
ground its velocity U relative to the inertial frame X is given by the Lorentz 
Transformation: [12 5] 

U = [(R + h) Q + V] / [1+ Q (R + h)V/c 2 ] (5.2.2). 

As in these flights c 2 » Rf2V [2] the we have U ~ (R + h)f2 + V; so the time 
interval dx 2 = dx A , measured by a fixed clock in the plane is given by 

(di A ) 2 = {[1-2GM/c 2 (R + h)] - [(R + h)Q + V] 2 /c 2 } dt 2 (5.2.3), 
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where the velocity V ~ 300 m/s has signals ± depending on the direction of 
movement of the airplane in relation to the ground. Measurements using 
high precision clocks have confirmed [2,7] the predictions given by (5.2.1) 
and (5.2.3). We recall that in SR the transverse Doppler effect was 
introduced "ad hoc" to explain the kinematic effect because the reference 
system where clock 1 is fixed was non-inertial with an radial acceleration 
a R = Q"R. From this we verified that this kinematic effect can only be 
explained rigorously within the context of the GR. 

(5.3) Gravitational & Longitudinal Kinematic Doppler Effect. 

Now let's consider that a light detector (clock) is fixed at a point P 0 
and that a light emitter at a point P is moving directly away or approaching 
P 0 (see Figure 1). The radial Cartesian velocity of P defined by q = ds/dt, 
according to (5.1), would be given by [3] 

q = {(dr/dt) 2 /(l-2y/r) + r 2 (d0/dt) 2 + r 2 sin 2 0(dcp/dt) 2 } l/2 /( 1 -2x/r) 1/2 = 

= (l-2%/r)‘ 1/2 {(l-2x/r) _1 (dr/dt) 2 + r 2 (d0/dt) 2 + r 2 sin 2 0(dcp/dt) 2 } 1/2 (5.3.1) 

Defining the radial coordinate velocity by V = dr/dt and assuming 
that d0 = dcp = 0 we get from (5.3.1): 

q = (l-2 X /r)- 1/2 V (5.3.2), 

showing that q —> V when r —>-go or M —> 0 in absence of matter. It can 
also be obtained [3] the general y relation , 

y 2 = (1 - 2%/r)/(l - q 2 /c 2 ) (5.3.3). 

We see that y = (l-V 2 /c 2 ) 172 = (l-(3 2 ) 172 , with (3 = V/c, when r — >oo or M —>0. 

Let us suppose that the from the point P are emitted two pulses of 



Figure l.Light detector fixed at P 0 and light emitter in motion at P. 
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light with own frequency v at the events E and E' defined in the coordinate 
system X by E = (t, r, 9, cp) and E'= (t + dt, r + dr, 0, (p). They are separated 
by the invariant proper-time interval 

ds 2 = (l-2%/r)c 2 dt 2 + dr 2 /( 1 -2%/r) 

= (l-2x/r)(l - q 2 /c 2 )dt 2 (5.3.4), 

where the radial velocity q is given by q = (dr/dt)/(l-2%/r), according (4.1.4) 
putting d(p/dt = d9/dt = 0. The reception of the first signal by (fixed) P 0 is 
defined by the event E 0 defined by E 0 = (t Q , r o ,0,0). The reception of the 
second signal is defined by E 0 ' = (t G + dt Q , r Q , 0, 0). These two events are 
separated by the interval ds G given by 

ds„ = (l-2x/r 0 ) 1/2 cdt 0 (5.3.5), 

and therefore the frequency, v 0 + dv 0 , and the wavelength, X 0 + dk 0 , 
received by detector P 0 are given by 

v 0 + dv G = c/ds 0 and X 0 + dk 0 = cds 0 (5.3.6). 

Since the events E & E 0 and E' & E 0 ' are connected by light signals they 
can be taken as null geodesics through P and P 0 putting ds = d9 = d(p = 0 in 
(5.3.1), obtaining (1-2% /r) c 2 dt 2 + dr 2 /(1-%/r) = 0 , that is, 

dr/dt = c(l-2%/r) (5.3.7). 

So, for the pair E&Eo we have, integrating (5.3.7) with t Q —> t and r () —> r : 

r 

t t 0 = (l/c) dr/( 1 -2 / 1 -) 

<b 

(5.3.8) . 

Similarly, for the pair E'& E a ' integrating (5.3.4) with t 0 -t-dt () —> t+dt and 
r 0 +dr 0 —» r+dr we get 

t + dt t 0 + dt 0 = (l/c) dr/(l-2 /r) 

'o+ dr o 

(5.3.9) . 
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Subtracting (5.3.8) and (5.3.9) results dt 0 - dt = - (1/c) dr/(l-2%/r) = -(q/c)dt, 
gives 

dt Q = (1- q/c)dt (5.3.10). 

In this way, from (5.3.5), (5.3.9) e (5.10) we obtain, 

d So /ds = {(1- 2 X /r„)/(l- 2 X /r) } m {(1 - q/c)/(l + q/c)} 1 ' 2 (5.3.11) 

Since v = 1/ds (5.3.11) can also be written as 

v/v„= {(1- 2 x /r 0 )/(l- 2 x /r)} 1/2 {(1 - q/c)/(l + q/c)} 1/2 (5.3.12). 

Note that v G is the frequency measured by a detector (P 0 ) fixed at the Earth 
and v is the proper frequency emitted by a system (P) in motion with the 
radial velocity q relatively to the Earth (P 0 ). 

From (5.3.12) we see that the gravitational and kinematic Doppler 
Effects are intertwined. (5.3.12) predicts only kinematic effects by making 
G = x = 0 and only gravitational effects when q = 0. 

Emitter and receiver fixed at points E and R. r81 

In this particular case the two points E = (t E ,r E ,9 E ,(p E ) and 
R = (t R ,r R ,0 R ,(p R ) are connected by a null geodesics defined by the equation 

(1-2%/r)c 2 1 +2 - (1-2%/r)" 1 r +2 - r 2 0 +2 - r 2 sin 2 0 (p +2 = 0 (5.3.13), 

where, generically, f - = df/du. From (5.3.13) we have obtain 

dt/du = (l/c){(l-2%/r)' 1 g ;j (dx‘/du)(dx J /du)} 1/2 , where = - gy . 

Integrating this expression we get 

t R - t E = (1/c) {{(1-2%/r)' 1 gy (dx 1 /du)(dx J /du)} 1/2 du. 

As this integral depends only on the path through space, so with a spatially 
fixed emitter and spatially fixed receiver, t R - t E is the same for all signals, 
land 2, sent. Thus, for two signals we have t R (1) - t E (1) = t R (2) - t E (2) 
showing that time intervals measured in E are equal to that measured in R, 
that is, 

AtR = t R ,2) - t t (1> = t E ,2) - t E n ’ = At E (5.3.14). 

That is, the coordinate time difference at the point of emission E equals the 
coordinate time difference at the point R of reception. However, the clocks 
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of an observer at R and of an emitter at E record proper times and not 
coordinate times. The record proper times At r and At e are given by [8] 

Ax R =(l-2x/r R ) 1,2 At R and Ax E =(l-2 X /r E ) 1,2 At E (5.3.15). 

Since At R = At E we have 

At r /At e = [(1 -2 X /r R )/( 1 -2 X /r E )] 1/2 (5.3.16), 

that is a particular case of (5.3.12) when q = 0, remembering that Av =1/Ax. 


(5.4) Particle Geodesics in Weak Gravitational Field Limit. 

Let us see recover the Newtonian equations of motion of a massive 
particle using as a starting point the geodesic equations in the limit of very 
weak gravitational field (O) and with velocities dxVdt « c (i =1,2,3). 
Suppose a system where the metric tensor is given by g^ v = r^ v + h^ v where 
r\x v = diagf-1,1,1,1] and h^is very small but cannot be neglected. In these 
way one can show, after a somewhat long calculation, that [8] 

md 2 x7dt 2 = - m8 IJ Oj(c 2 hoo/2) + mc5 lk (<3jh ok - 8 k h oj )(dx j /dt) (5.4.1). 

The first term on the right would be the gravitational force -m grad(O) 
where ® = c h 00 /2 + const and the second term would be a Coriolis force 
that is null in a nonrotating reference system. Taking const = 1 we see that 
goo in the Newtonian approximation would be given by g OQ = 1 + 2®/c 2 . 

Example 1.Newton's Law of Universal Gravitation. 

Using (5.4.1) we recover Newton's law remembering that g 00 in the 
SM in the weak field limit (for large l/%) in given by g 00 = (1-2% /r) ~ 1 - 
2GM/rc 2 , that is, h 0G = - 2GM/rc 2 . This would imply that V = - GM/r and 
consequently, using (5.4.1): 

F(r) = md 2 r/dt 2 = -m grad(V) = - (GMm/r 2 ) r 

in agreement with the gravitational Newton's law. 

(6) Principle of Equivalence and Geodesic Principle. 

In many basic textbooks [1 ' 3,5] about the Theory of Relativity are 
found analysis on the Principle of Equivalence. This is proposed as an 
immediate consequence of recent experimental measurements [2] which 
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show that the gravitational mass m g is equal to the inertial mass m, with an 
accuracy of one part in 10 9 . We suggest the lecture, for instance, of the 
analysis presented by Yilmaz [5] where is shown Figure 2: 



Figure 2. A system (a) which is stationary in a uniform gravitational field g is 
physically equivalent to a system which is in a gravitational-free space (b) but 
accelerated in the opposite direction with an acceleration y = g. 

The Principle of Equivalence can be written, for instance, as: 

(6.a) "It is impossible to tell whether a system (room) is in a state of 
accelerated motion by experiments performed in that system alone. For the 
observer in the room can claim that this room is stationary but there is a 
gravitational field present in the room." 

(6.b)"A system which is stationary in a gravitational field of strength g is 
physically equivalent to a system which is in a gravitation-free space but 
accelerated in the opposite direction with an acceleration g.” 

(6.c)" In physical terms the effects of a gravitational field can be removed, 
locally at least, by employing an appropriately chosen accelerated 
coordinate system. This is a local coordinate system which is falling freely 
in the gravitational field of the distribution of matter". 

According to Section 4 it is always possible to determine a Locally 
Inertial coordinate-system (LIC) at any point of a Riemannian space. The 
metric at the origin of this coordinate-system is of aflat space (Appendix A). 
This implies that the Riemann-Christoffel R\ vli , the Ricci R > 4 , and the 
Curvature Invariant R tensors are equal to zero (see Appendix A). In any 
flat Riemannian space-time (x), representing a large distribution of matter, 
the motion of a small particle will be described by the time-like geodesics 
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d 2 x7ds 2 = 0 (o=l,2,3,4) (6.1), 

which represent the motion of a particle locally identical with the motion 
under Newton's Fist law. They are therefore immediately integrable to give 
dx°/ds = v° and x°- x 0 ° = s v° where the four v° and the four x 0 ° are 
constants of integration. In physical terms, these LIC may be said to form a 
system which is falling freely in the gravitational field of the distribution of 
matter which corresponds to the formulation (6.c) of the Principle of 
Equivalence. In other words, the existence of a LICS at every event in 
space-time expresses the Principle of Equivalence and indicates that the 
effect of a gravitational field can be abolished locally by the choice of a 
suitable accelerated coordinate-system. The representation of the history of 
the motions of small test-particles by time-like geodesics and of light-rays 
by null-geodesics (see Section 3.3) in Riemannian spaces form the 
"Geodesic Principle". As stated above, the geodesic principle presupposes 
that the distribution of matter is given through its energy-tensor, that the 
appropriate Riemannian space-time has been determined through Einstein's 
equations and that the particle itself contributes nothing to the distribution 
under whose gravitational influence it moves. [3] See comments of 
McVittie [3] about geodesic principle, continuous distributions of matter, 
cosmology and Einstein's equations. 

(7)Accelerated Coordinate-Systems in Special Relativity. 

It has been mentioned above that the existence of a LICS at every 
event in space-time expresses the Principle of Equivalence and indicates 
that the effect of a gravitational field can be abolished locally by the choice 
of a suitable accelerated coordinate-system. This result has led to some 
confusion for it has been interpreted as equivalent to the statement that the 
essence of GR consists in the use of accelerated coordinate-systems, in 
contrast with SR, where relatively non-accelerated coordinate-systems are 
employed. [3] This statement is, however, misleading because there is 
nothing to prevent the investigator from using accelerated coordinate- 
systems is SR if he chooses to do so. Let us consider, for instance, an 
inertial system S in the Minkowski space-time (cT,X,Y,Z), with metric 

ds 2 = c 2 dT 2 - dX 2 - dY 2 - dZ 2 (7.1). 
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In S the motion of a non accelerated the test particle would be governed by 
the time-like geodesics equations d 2 X^/ds 2 + I\/ (dXVds) (dX v /ds) = 0, 
where the Christoffel symbols = 0, that is, 

d 2 X^/ds 2 = 0 (p=l,2,3,4). (7.2). 

The Riemann-Christoffel and the Ricci tensors in S are identically zero and 
so is the energy-tensor. 

Consider now a non-inertial Minkowski space-time coordinate- 
system S' in which the (ct,x,y,z) are obtained by Lorentz transformations 
from (cT,X,Y,Z). These transformations do not alter the null-character of 
the Riemann-Christoffel and the Ricci tensors of the MST which remains a 
flat space-time (see Appendix A). Hence the Ricci tensor is still identically 
zero and so is the energy-tensor, in spite of the fact that the coefficients 
g ik (xj) of the new metric -ds 2 = g ik (xO dxi dx k can be no longer constants. 
Thus, ds 2 does not represent any distribution of matter and therefore no 
gravitation field is present. Note that the equations of motion of a particle 
in S' can be difficult to obtain in, general case, when the Lorentz y ^ 1. 

An example of this can be seen, for instance, in reference [3], pag.78. In 
Sections (7.1) and (7.2) are shown simple cases when y = 1. 

(7.1) Rotating Reference System. 

(7.1.a) Cartesian Coordinates. 

Let us pass from an inertial reference frame S = (T, X, Y, Z) to a 
non-inertial rotating referential system S' with constant angular velocity co 
along the z-axis. 



Figure 2. Coordinate system S' = (ct, x, y, z) rotating relative to the coordinate system 
S = (cT, X, Y, Z). 

Defining a new referential system S' = (t, x, y, z): 


22 



T = t, 

X = x cos(cot) - y sin(cot), 

Y = x sin(cot) + y cos(cot), 

Z = z, 

we see that ds 2 = c 2 dT 2 - dX 2 - dY 2 - dZ 2 now becomes written as 

ds 2 = c 2 dx 2 = [c 2 - co 2 (x 2 +y 2 )]dt 2 + 2cox dy dt - dx 2 - dy 2 - dz 2 (7.1.al). 

For S' the geodesic equations are d 2 x T /di 2 + r>. V T (dxVdx) (dx v /dx) = 0, 
where r >A , T = (l/2){ Sg-J6x v + dg Vil /dx x - og 4 ,/fx v }, (x = 0,1,2,3). Since 
goo = [ c 2 - co 2 (x 2 +y 2 )], g o2 = 2cox and g u = g 22 = g 3 3 = - 1 we see that the 
only non null r^ v T are r 0 o = - co x, r 20 = -2co, r 0 o“= - co y and T i 0 = 2co. 
Thus, the equations of motion for a free massive particle are given by, [8] 

d 2 t/dx 2 = 0 , 

d 2 x/dx 2 - co 2 x(dt/dx) 2 - 2co(dy/dx)(dt/dx) = 0, (7.1 .a2) 

d 2 y/dx 2 - co 2 y(dt/dx) 2 + 2co(dx/dx)(dt/dx) = 0, 
d 2 z/dx 2 = 0. 

As the first equation implies that dt/dx = k = constant, that is, dt = k dx the 
remaining (7.1.2) equations may be written as 

d 2 x/dt 2 - co 2 x - 2co(dy/dt) = 0, 

d 2 y/dt 2 - co 2 y + 2co(dx/dt) = 0, (7.1.a3), 

d 2 z/dt 2 = 0 

These calculations clearly shown that in the non inertial system S' the 

acceleration has a geometrical nature. Multiplying by m the noninertial 
accelerations given by Eqs.(7.1.a3), written in 3-dim vector notation, we 
get the noninertial forces 

md 2 r/dt 2 = m co x (co x r) - 2mco x (dr/dt) (7.1.a4), 

where r = (x,y,z) and co = (0,0,co). The component 2mco x (dr/dt) is the 
centrifugal force and - 2mco x (dr/dt) the Coriolis force. [1,9] . These forces 
are not due to physical interactions. Note that with the choice T = t, 
where T is the proper time measured by clocks at rest in S, the time t is 
exactly the proper time for an observer situated at the common origin O 
of the two systems. So, observers close to O who are at rest in the rotating 
system would accept Eqs.(7.1.a3) and (7.1.a4) as approximately valid and 
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recognizes the terms on the right as forces of acceleration.^ 1 

Similar results can be deduced using the "weak field approximation" 
seen in Section (5.4), when we have assumed that g^ v = + h^ v where 

r^v = diagf-1,1,1,1] and h^is very small but cannot be neglected. The terms 
—► - co 2 (x 2 +y 2 )/c 2 , coy/c, - cox/c would be small for small angular 
velocities co and so, (7.1.a3) and (7.1.a4) can be obtained. [8] 

(7.1.b)Cylindrical Coordinates. 

Transforming Cartesian (t,x,y,z) to cylindrical coordinates putting 
z'= z, x = p cosq), y = p sincp from (7.1 .a2) we get the "accelerations" 

d 2 t/dx 2 = 0 , 

d 2 p/dx 2 - co 2 p(dt/di) 2 - p(dcp/dt) 2 - 2cop(d(p/dx)(dt/di) = 0, 

d 2 (p/di 2 - co 2 p' 1 (dp/dT)(d(p/dT) + 2ojp _l (dp/di)(dt/di) = 0, (7.1.bl) 

d 2 z/di 2 = 0. 

Introducing the mass m of the particle and rearranging we have the virtual 
forces, noting that dt = k di, 


(A) m[d 2 p/di 9 ] = p[(d(p/dt) 2 + mco 2 ] + 2mcop(dcp/di), 

(B) m[d 2 (p/di 2 ] = - 2m(dp/dT)(dcp/di)/p - 2mco(dp/di)/p, (7.1.b2). 

(C) m(d 2 z/di 2 ) =0, 

The first term of (A) is the centrifugal force due to particle's own motion 
and of the rotation of the frame and the second term is the radial Coriolis 
force. The right hand side of (B) is the tangential Coriolis force; it depends 
on the radial velocity, and on the angular velocity co of the particle in the 
frame and rotation of the frame. 

Now, if the inertial cylindrical coordinate system is represented by 
S'= (t, p',cp', z') we have 

ds 2 = c 2 dt 2 - dp' 2 - p' 2 dcp' 2 - dz' 2 (7.1.b3). 

If t, p, cp and z are the cylindrical coordinates of the rotating system 
S = (t, p,cp, z) we see that p'= p, z'= z and cp'= cp + cot. In this way (7.1.b3) 
becomes given by 

ds 2 = (c 2 - co 2 p 2 )dt 2 - 2cop 2 dcpdt - dz 2 - p 2 dcp 2 - dp 2 (7.1.b4). 
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(8) Distances and Time Intervals in GR. 

In GR the choice of a reference system is not limited by nothing; the 
3 coordinates x',x 2 , x 3 can be arbitrary quantities defining the position of 
the particles in the space and the temporal coordinate x° can be determined 
by a clock marking its proper time. The main problem is how one can 
determine using x^x^x 2 , x 3 the real distances f and the real time lapses x. 

(8.a) Real time intervals. 

Let us first determine the connection between real times x and the 
coordinates x°. To do this let us estimate the interval ds between two 
events that occur in the same point, that is, when dx 1 = dx 2 = dx 3 = 0. So, 
as - ds 2 = g ik dxidx k we get ds 2 = - c 2 dx 2 = g 0G dx G 2 , that is 

dT = (l/c)(-g 0o ) 1,2 dx" (8.a.l), 

from which we see that the finite real time interval x between two arbitrary 
events occurring in the same point of the space is given by, 

T = (l/c)^(-g 00 ) 1,2 dx° 

(8.a.2). 

The relations (8.a.l) and (8.a.2) show how to determine the real 
times x (that is, the proper times in one point of the space) as a function of 
the coordinate x G . Note that, to have physical meaning, we must have 

goo < o (8.a.3). 

It is important to stress that the spatial components of tensor g,j must be 
positive and that g 00 must be negative. One metric tensor that does not 
satisfy the last condition (8.a.3) cannot represent a real gravitational field. 
If (8.a.3) is not obeyed would signify that the correspondent reference 
systems could not describe real bodies. [1] 

Weak gravitational field. 

As will be seen in Section (5.4) in weak gravitational fields ®(r) we 
have at a given point of the space, putting x°/c = t: 

X = (x°/c)(- g G0 ) 1/2 ~ (x%) [1+ 2 ®/c 2 ] 1/2 ~ (1 + ®/c 2 ) t (8.a.4). 
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When <J> < 0 the real or proper time x flows more slowly than the absolute, 
universal or coordinate time t = (x°/c) which is measured in absence of 
gravitational fields. 

In the Newtonian approximation, according to Section 5, g QO in 
would be given by g G0 = (1-2%/r) = 1 + 2®/c 2 ~ 1 - 2GM/rc 2 . If, for one is 
measuring the time fixed on the Earth surface, taking G ~ 6,7 10' 11 MKS, 

M ~ 6 10 24 kg, r ~ 6,4 10 6 m and c ~3 10 8 m/s using (8.a.4) he will see that 
x/t ~ 1-10 . That is, the time measured on the earth flows more slowly 
than the time measured in absence of the gravitational field. 

Time Synchronization. 

Now let us present the notion of Simultaneity in GR. That is, let us 
see if it is possible to synchronize watches that are in two different points 
of the space, in other words, to obtain a correspondence between their 
indications. Let us consider the general case when the metric tensor 
depends on the time coordinate x°. Such synchronization must be done, 
evidently, by an exchange of luminous signals between these two points. 
To do this let us consider the process of propagation of light signals 
between two infinitelly close points A and B represented by Figure 3. If x Q 
is the time of arrival of the signal in A, the moments of its departure from B 
and its arrival in B will be, respectively, x Q + dx Q (1) and x Q + dx 0 (2) . 

x °+dx a (2) 

x°+dx a f,) 

4 B 

Figure 3.Propagation of ligth signals between two infinitesimally close points A and B. 

In Figure 3 the parallel straight lines represent the universe lines 
corresponding to the coordinates x“ and x a + dx a , and the dotted sections 
represent the universe lines of the luminous signals(to simplify we can 
suppose that dx G (1) < 0 and dx 0 (2) > 0). Putting - ds 2 = g a p dx a dx p + 2g oa 
g 00 dx°dx a + g 0 odx 0 2 we have 



dx 0 (1) = -{ g oa dx a - [(goa gop - gapgoo) dx“dx P ] 1/2 }/g 00 
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dx 0 (2) = -{ goadx“ +[(g oa gop - gapgoo) dx a dx p ] 1/2 }/g GO (8.a.5) 

It is clear that the total lapse of "time" between the emission and the return 
of the signal to the same point B is equal to 

dx 0 (2) - dx 0 (1 ’ = -(2/goo) [(goa go(3 - gapgoo) dx a dx p ] 1/2 ( 8.a.6). 

Let us stipulate as simultaneous with the instant x° at the point A the 
indication x° of the clock at B located at the half of the instants of emission 
and return of the signal in this point, that is 

X°(B) = x° + Ax° = x°(A) + (l/2)[dx 0 (2) + dx 0 (1) ] (8.a.7). 

In other words, %°(B) is "simultaneous" with x°(A). Using (8.a.6) we obtain 

Ax° = - g oa dx a /g 00 (8.a.8), 

which is the "time" difference between two simultaneous events happening 
at points infinitesimally distant A and B. Since t = cx° (8.a.7) can also be 
written, according to Einstein, as 

t(B)= t, A +(t 2 A -t, A )/2 = (t 2 A + t, A )/2 (8.a.9), 

where t 2 A = c(x° + dx 0 ' 2 ') and t] A = c(x° - dx 0 (l) ). 

In stationary fields that is, when the metric tensor depends on the 
coordinate time x° the clocks synchronization in all space is impossible. 
For static fields, that is, when the metric is time independent this is 
possible. In this case, using (8.a.8), we can determine the simultaneity of 
two events performing the integration 

Ax °=- $go« dxa/ goo 

(8.a.l0). 

In stationary fields is possible to synchronize clocks only along a closed 
path, returning to the initial starting point, obtaining the time difference At 
by the cyclic integration^ ] of (8.a. 10): 

At = - (l/c)f g 0o dx" / g 00 

(8.a.ll). 
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In particular, in the case of a rotating coordinate system, taking into 
account the metric given by (7.1.b4), we have, 

At = (l/c 2 )J cop 2 d(p /(1 -co 2 p 2 /c 2 ) (8.a.l2). 

(8.b) Real Spatial distances. 

In SR one can define the element d£ of real spatial distance as the 
interval between two events spatially close dx° (o =1,2,3) happening at the 
same "instant" x°, that is, putting dx° = 0 in - ds 2 = g?, v dx\ix v . In GR it 
cannot in general be done because the proper time in a gravitational field is 
differently connected with the coordinate x° in different points of the 
space, according to (8.a.l) and (8.a.2). The calculation of the element df 
is given, for instance, by [U ] 

df 2 = y a (3dx a dx p and y a(3 = gap - (g 0 « g 0 p)/goo (8.b.l), 

where in the product dx a dx p we have only a and (3 = 1,2,3. That is, 
temporal components do not contribute. 

Example. Cylindrical Coordinates Rotation. 

Using the cylindrical metric shown in Section 7 and (8.b.l) we get 

ds 2 = -(1 - X 2 ) dx o 2 + dxi 2 + (2px)dx 2 dx G + x^dxo 2 + dx 3 2 (8.b.2), 

where ct = x 0 , p = x 1? cp = x 2 e z = x 3 and % = cop/c. Since the product 
dx 2 dx 0 are not taken into account we have the non null terms: 

goo = -(1 - X 2 ), gll= grr = 1, g22 = g«pq> = p 2 and g 33 = g zz = 1. 

In this way, from y a(3 = g a(3 - (g oa g 0 ( 3 )/g 0 o we get yn = gn-(goi goi)/goo= 1, 
y 22 = g22-(go2 go2)/goo= p 2 /(l - X) and y 33 = g 33 -(go3 go 3 )/goo = 1. So, 

df 2 = yii dx'dx 1 + y 22 dx 2 dx 2 + y 33 dx 3 dx 3 = dp 2 + dz 2 + y 22 d(p 2 . 

Finally, we can write 

df 2 = dp 2 + dz 2 + p 2 dcp 2 /(l - co 2 p 2 /c 2 ) (8.b.3). 

From (Bb.3) we see that the "real" circumference of a circle with radius p is 

f = 27tp/(l - co 2 p 2 /c 2 ) 1/2 > 27ip (8.b.4), 
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where 2np is the circumference measured by one that is not in the rotating 
system. So, if one fixed on the Earth is measuring the circumference of the 
equator, putting p ~ 6.300 km, c ~ 300.000 km/s , co = 2 tt/T and T = 24 h ~ 
9 10 4 s will verify that f/27ip = 1/(1 - co 2 p 2 /c 2 ) 1/2 ~ 1+10' 13 . 


Appendix A. Curvature Tensor and Space Curvature. 

In GR there are two important tensors, one of rank four Recalled 
Riemann-Christojfel tensor and another of rank two R > 4 , named Ricci 
tensor. The first one is defined by 

RVn - d(T^)/dx» - 0(V)/3x v + (rO(V) - (rV)(iV) (A.i) 

and the second one obtained by index contraction of (A.l): 

R ^ = d(T Xa ydx» - a(V)/3x° + (r, 0 T x V) - (iyxr IO °) (a.2). 

The Riemann-Christojfel tensor (A.l) is known as the curvature tensor and 
the scalar 

R u = g te R„x (A. 3) 

is named curvature invariant. It measures a property of the Riemannian 
space that it is analogous to the curvature of a 2-dim surface described by 
(2.1.5) ds 2 = a 2 (d0 2 + sin 2 0 d(p 2 ). Putting 0 = x 1 and cp = x 2 we verify that [3,10] 

R u = _ 2/a 2 (A.4). 

As the Gaussian curvature [3,10] of a sphere of radius a is known to be 1/a 2 
we verify that (A.4) shows that curvature invariant R ' is proportional to 
the Gaussian curvature. Note that when the space dimension is larger than 
2-dim the relationship of these tensors to the "curvature" of the space is 
even more remote. 

In certain Riemannian spaces of n-dim it is possible to find a 
coordinate-system (X) covering all the points of the space, in terms of 
which the metric becomes 

ds 2 = I, £l (dX 2 ) 2 (A.5) 

where the & x are positive or negative constants. In these conditions one can 
show [3] that (A.l), (A.2) and (A.3) are equal to zero. Such spaces ar eflat by 
analogy with the Euclidean space whose metric is of the type (A.5) when 
Cartesian coordinates are used. 

In a flat space one has necessarily R\ V)i = 0. It is possible to show the 
converse; if R\ V)i = 0 throughout the space, then the space is flat. In other 
words, the space is flat if the curvature tensor R\ V)J is zero. 
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Appendix B. Absolute Differentiation, Covariant Derivative 

and Geodesics. 

(B.l) Absolute Differentiation. 

The absolute differentiation^' 3,5,91,1 ^ of a vector field A*\ for instance, 
is denoted by the symbol D defined by the operation 

DA 11 = dA 11 + IV A x dx v (B. 1.1), 

where d is the ordinary differentiation. 

(B.2) Covariant Derivative. 

The covariant derivative along a curve is defined by 

DATDx = A p ; v (dxVdi) = dA7dx + Yff A L (dx v /dx) (B.2.1), 

which is a vector because it is a product of the tensor A* 1 . v and the vector 
(dx v /di). Note that dA7dx is not a vector. [2] 

Along a geodesic the velocity field v 11 = dx7dx which is tangent to 
the curve in each point obeys the equation 

dV'/dx 2 + T> V T (dxVdx) (dx v /dx) = 0, that is, 

dv7dx + WV = 0 

This implies, using (B.2.1), that 

Dv7dx = 0 (B.2.2), 

that is, one verifies that along a geodesic v^ = constant. Note that in a 
Riemannian space ("curve") as W 0 a particle even in absence of a 
physical interaction has inertial accelerations^^ Section 7). 

For flat spaces found in Special Relativity(SR) a free particle, not 
submitted to external forces, moves along a geodesic ("straight line") with 
constant velocity v* 1 and null acceleration dVVdt 2 . As a null geodesic in 
Minkowski Space (MS) is given by ds 2 = (dx 4 ) 2 - (dr) 2 = 0 we get (dr/dt) 2 
= v 2 = c 2 , showing that it represents the history of a light ray motion. 

In GR in an obvious generalization of the SR, it is postulated that a 
free particle describes a geodesic. The hypothesis that free particles 
describe geodesics is called Geodesic Principle and corresponds to an 
extension of the Inertia Principle of Galileo. In fact in GR the "free 
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particle" is immersed in a Riemannian geometry dictated by a gravitational 
field. The interaction is "substituted" by the geometry. 

Appendix C. Spacetime Intervals and Vectors. 

In Euclidean space the separation between two points is measured 
by the distance between two points. This distance is purely spatial and is 
always positive. In a Riemannian spacetime the separation between two 
events takes into account not only the spatial separation between the events 
but also their temporal separation. In this case, between two events, we 
have spacelike intervals that are defined taking into account the metric 

ds 2 = g^ v dxMx v (C.l). 

(a) ds 2 < 0 —> Timelike interval. 

2 2 2 

In Minkowski spacetime we have c At“ > Ar . For two events 
separated by a timelike interval there is enough time between them to exist 
a cause-effect relationship between them. 

(b) ds 2 > 0 —> Spacelike interval. 

In Minkowski spacetime we have c 2 At 2 < Ar 2 . When a spacelike 
interval separates two events there is no enough time between their 
occurrences that can be explained as being created by a causal relationship 
crossing the spatial distance between the two events with the speed of light 
or slower. 

(c) ds = 0 —> Lightlike interval or null interval. 

In Minkowski spacetime c 2 At 2 = Ar 2 . In lightlike interval the spatial 
distance between two events is exactly the time interval elapsed between 
the two events. Events connected by photons have lightlike or null 
intervals. 

Vectors. 

In Riemannian spacetime a vector q A is defined as timelike, spacelike 
or null if the following conditions are, respectively, obeyed [8] 

gx, v (x) qV < 0 , > 0 or = 0. 
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Appendix D. Paths of particles moving in the vicinity of spherical 
massive objects using Lagrangian formalism. 

As seen in Section (3.4) the path of massive particles and photons 
described, respectively, by timelike and null geodesics can be obtained 
using the Lagragian formalism. This is done solving the Euler-Lagrange 
equations 

d(5L/ox M )/ds - oL/ox M = 0 (D. 1), 

where x a = x a (u), x a = dx a /ds and 

L(x a , x a ) = (l/2)g^ v (x a ) (dx7ds) (dx v /ds) (D.2). 

Solving (D.l) let us see how to determine the motion of massive 
particles and of photons in the vicinity of a spherical massive object. 
Following Sec. 5 the metric for a spherical object with mass M is the SM 
ds that in spherical polar coordinates (r, 0, cp) is given by 

ds 2 = (l-2x /r) c 2 dt 2 - dr7 (l-2x/r) - r 2 d0 2 - r 2 sin 2 0d 2 (p (D.3), 
where % = GM/c 2 , x° = t, x 1 = r, x 2 = 0 and x 3 = (p. 

(D.a) Massive Particle. 

Taking into account the metric (D.3) one can show that the 
Lagrangian (D.2) becomes written as [8] 

L(x“, x“) = (1/2){ c 2 (l-2x/r)(dt/ds) 2 - (l-2x/r) _1 (dr/ds) 2 - r 2 [(d0/ds) 2 + (dcp/ds) 2 sin 2 0 } (D.a.l). 

Because of the spherical symmetry, there is no loss of generality in 
confining our attention to particles moving in the "equatorial plane" with 
0 = 7i/2. In this case (D.a. 1) is simplified and for p = 2, that is, x = 0 we 
get from (D.l) and (D.a.l), 

F(r)"'(d 2 r/ds 2 ) + (mc 2 /r 2 )(dt/ds) 2 - F(r)' 2 (m/r 2 )(dr/ds) 2 - r(dcp/ds) 2 = 0, (D.a.2) 
where F(r) = (l-2x,/r). 

Since t and (p are cyclic coordinates we get from (D.l), putting 
t + = dt/dx and (p + = dcp/dx : 

5L/5t + = const and 5L/5cp + = const (D.a.3). 

With (D.a.3) and 0 = tt/ 2 in (D.a.2) we obtain 

(l-2x/r)t + = k and r 2 (p + = h (D.a.4), 

where k and h are integration constants. The first one gives the relation 
between the coordinate time t and the proper time x and the second one is 
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analogous to the equation of angular motion conservation. We have from 
(D.3) with 9 = nl2\ 

c 2 F(r)t +2 - F(r) V 2 - r 2 (p +2 = c 2 (D.a.5), 

which is much less complicated than (D.5) and yields an equation 
analogous that expressing the energy conservation. Indeed, taking into 
account (D.a.4) we obtain, putting u =l/r and where E =c 2 (k 2 -l)/h 2 : [8] 

(du/dcp) 2 + u 2 = E + (2MG/h 2 )u + (2GM/c 2 )u 3 (D.a.6). 

The analogous Newtonian equation which is given by 

(du/dcp) 2 + u 2 = E + (2MG/h 2 )u (D.a.7), 

where h = r 2 (dcp/dt) is the angular momentum per unit of mass has a well 
known solution [8,9] 

u = 1/r = (GM/h 2 )[l + e cos(cp - cp 0 )]. 

The general relativistic effect introduced by the cubic term u in 
planetary motion in (D.a.6) is analyzed in reference 8 (pag.144-146). 

Vertical Free-Fall. 

In this case cp = constant, thus from (D.a.6) —> cp + = 0 —> h = 0. In this 
way putting cp + = 0 and (1-2% /r) t + = k in (D.a.5) we obtain, 

r +2 - ck 2 +c 2 (l-2x/r) = 0 (D.a.8). 

If the particle is at rest (r + = 0) when r = r 0 —► k 2 = 1- 2x/r 0 , showing that k 
is not an universal constant, but depends on r 0 , that is, on the geodesic. In 
particular, if r + —> 0 when r —> oo , then k = 1. 

Differentiating (D.a.8) we get, 

r ++ + GM/r 2 =0 (D.a.9). 

This equation has the same form ( d 2 r/dt 2 + GM/r 2 = 0) as that found in the 
Newtonian approach; remembering that in (D.a.9) the coordinate r is not 
the vertical distance and the derivative r ++ is with respect of the proper 
time x of the particle in motion, not the universal time t. 

Putting k 2 = 1- 2%/r 0 in (D.a.8) we obtain 

r +2 /2 = MG(l/r - l/r G ) (D.a.10) 

which is "similar" to that in the Newtonian theory. Integrating (D.a.10) 
assuming that r G is the initial coordinate of particle at x = 0 we can calculate 
the proper time x experienced by the particle in falling from r 0 up to r : 

= (1/V2GM) ° [r 0 r/(r 0 r)] 1/2 dr (D.a.ll). 

r 

It permit us to establish the relation r = r(x). To obtain r = r(t) we take into 
account dt/dr = (dt/dx)/(dx/dt) with (D.a.4), 
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dt/dx = k/(l - 2yjr) = (\-2ylv 0 ) m l(l - 2yjr). 
and with (D.a.10) getting, with % = GM/c 2 : 

v(t) = dr/dt = -[r 3/2 (r 0 - 2m) 1/2 ]/{c(2m) 1/2 (r - 2m)(r 0 - r) 1/2 } (D.a.12) 

The way in which the coordinate time t depends on r for a radially 
falling particle becomes more comprehensible if, considering for simplicity 
that it is at rest at infinity r 0 —>oo, we compare its coordinate speed, 
obtained from (D.a.12), 

v(r) = Idr/dtl = (2mc 2 ) 1/2 (r - 2m)/r 3/2 

with the classical Newtonian speed [9] 

v c (r) = V2mc 2 /r 1/2 . 

This difference is clearly seen in Figure (D.l). 



Fig.(D.l). Comparison of the coordinate speed v(r) with the classical Newtonian speed v c (r) for 
a particle falling from rest at infinite. 

Motion in a circled 81 

For a circular motion in the equatorial plane we have r = const and, 
consequently, r + = r ++ = 0. Equation (D.a.2) then reduces to 

mc 2 t +2 = r 3 (p +2 giving (dcp/dt) 2 = GM/r 3 (D.a.13) 

From (D.a.13) we see that for Acp = 2 tt the coordinate time At for a 
complete revolution is given by 

At = 27i(r 3 /GM) 1/2 (D.a.14). 

This expression is exactly the same as the Newtonian expression for the 
period of a circular orbit of radius r, that is, Kepler's third law. However, in 
the relativist case we cannot say that r is the radius of the orbit, but we see 
that the spatial distance traveled in one complete revolution is 27ir, just as 
in the Newtonian case. 
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Observer 
alr = r 0 


Fig.(D.l). Spacetime diagram shown a circular orbit as viewef by a fixed observer at r = r Q . 

Fig.(D.l) is a space time diagram illustrating one complete revolution as 
viewed by an observer fixed at a point r = r 0 . B i is the event observed when 
the particle is at the start point of orbit at Ai, while B 2 is that of his viewing 
its completion at A 2 . If the coordinate time difference between Ai and A 2 is 
At the coordinate time difference between B1 and B2 is also At, according 
to (5.3.14). So, the proper time Ax 0 which the observer measures for the 
orbital period At is given by 

AT o =(l-2x/r 0 ) 1,2 At (D.a.15). 

As r 0 —>oo, Ax 0 —>oo, so At is the orbital period as measured by an observer 
at infinity. So, as At is directly observed, according to (D.a. 14) the 
coordinate r also could be measured if the mass M is known. 

Now, let us determine the proper lapse of time measured by an 
observer moving with the rotating particle A. In this case the relationship 
between t and x is obtained using (D.a.4) and (D.a. 13). So, 

t +2 = kV/(r - 2yf [a] and 

cp +2 = mc 2 t +2 /r 3 = mc 2 k 2 /[r (r-2 j) 2 ] [b] (D.a. 16). 

Substituting (D.a. 16) [a] in (D.a.5) and putting r + = 0, since now the 
observer moves in a circle, we obtain 

k 2 = (r-2 X ) 2 /[r(r - 3 X )] (D.a.18). 

Putting this k 2 in (1-2%/r) t + = k, given by (D.a.4), we get, as t + = dt/dx: 

Ax = (1-2%/r) k _1 At = [(r-3 X )/r] 1/2 At = 2jr{(r 3 /MG) (l-3MG/rc 2 )} 1/2 (D.a. 19). 

This equation shows that we can have circular orbits only for k 2 > 0, 
that is, only when r > 3% = 3GM/c 2 . In the limit r —> 3%, Ax —► 0, 
suggesting that photons can orbit at r = 3% .This indeed occurs(see reference 181 ). 
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Appendix E. Length and Time in Schwarzschild Spacetime(SST). 

The SST has the line element given by, with % = 2M/c 2 : 

ds 2 = c 2 dx 2 = (1-2% /r) c 2 dt 2 - dr 2 / (1-2%/r) - r 2 d0 2 - r 2 sin 2 0d 2 (p (E.l). 

If we take a slice given by t = constant we obtain a 3-dim manifold. Putting 
the line element 

ds 2 = g* ij dx 1 dx J (i,j = 1,2,3, x l = r, x 2 = 0, x 3 = cp ), g*y = - g,, (E.2), 

describing a slice which is a space rather than a spacetime. As g*y is time 
independent we can observe two events at the same point in space 
occurring at different times. When M = 0 the line element (E.l) describes a 
flat spacetime, while the line element (E.2) describes the Euclidean flat 
space in spherical coordinates according to (2.1.4). 

What about radial distances given by 0 and cp constants? The line 
element (E.l) shows that for these infinitesimal radial distance dR is given 
by dR = dr/( 1-2%/r) 1/2 , so dR > dr and r no longer measures radial 
distances. Not even the origin of the r coordinates can be defined! 
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